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A  FTiel  Engine  Evaluation  System  (FEES)  has  been  designed  to 
evaluate  the  effects  of  changes  in  fuel  composition  on  engine 
performance  and  operability.  This  report  describes  FEES  and  the 
tests  carried  out  to  calibrate  the  system  and  to  determine  its 
reliability.  From  the  results  obtained  recommendations  are  made 
to  improve  the  system.  y 

FEES  was  designed  to  handle  spark  ignition  and  compression 
ignition  research  engines  of  power  output  up  to  22  kW  at  3500  RPM. 
Cold  start  capability  testing  to  -40°C  is  also  available.  The 
engine  used  for  the  above  tests  was  a  Sachs,  single  rotor  KM914A 
rotary  combustion,  spark  ignition  engine. 


RESUME 


Un  systeme  d' evaluation  de  moteur  en  fonction  du  carburant 
(Fuel  Engine  Evaluation  System)  (FEES)  a  ete  congu  pour  evaluer  les 
effets  causes  par  des  modifications  de  la  composition  du  carburant 
sur  le  fonctionnement  et  la  performance  d'un  moteur.  Le  rapport  donne 
une  description  du  FEES  et  des  essais  d'etalonnage  et  de  fiabilite  du 
systeme.  A  partir  des  resultats  obtenus,  des  recommandations  sont 
formulees  pour  ameliorer  le  systeme. 

Le  FEES  a  ete  congu  pour  des  moteurs  de  recherche  a  allumage 
par  etincelle  et  a  allumage  par  compression  d'une  puissance  d'au  plus 
22  kW  a  3500  tr/min.  II  comprend  aussi  des  essais  de  demarrage  a  froid 
realises  jusqu'a  -40°C.  Le  moteur  utilise  dans  les  essais  mentionnes 
precedemment  etait  un  moteur  rotatif  Sachs  KM914A  monorotor,  a  allumage 
par  etincelle. 
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1 . 0  INTRODUCTION 

A  subprogram  to  investigate  Fuel/Powerplant  interaction  has  been 
initiated  at  the  Defence  Research  Establishment  Ottawa  (DREO) .  One  of  the 
central  objectives  of  the  subprogram  is  to  determine  the  effects  of  the 
anticipated  degradation  in  quality  of  future  domestic  and  foreign  fuels  on 
Canadian  Forces  powerplants.  An  intramural  study  has  been  implemented  to 
complement  extramural  activity  for  the  investigation  of  compression  ignition 
powerplants  combusting  future  Canadian  fuels. 

To  monitor  engine  performance  and  operability  while  combusting 
fuels  with  low  cetane  number  and/or  wider  boiling  range,  and  fuels  derived 
from  tar  sands  blends,  a  Fuel  Engine  Evaluation  System  (FEES)  has  been 
designed  and  assembled.  This  report  endeavours  to  discuss  this  system  and 
the  tests  performed  during  its  shakedown  and  preliminary  calibration  trials 
using  the  KM91'*A  Sachs  Rotary  Combustion  Engine  as  a  prime  mover. 

2 . 0  TEST  OBJECTIVES 


The  tests  described  in  this  report  were  conducted  to  fulfill  the 
following  objectives: 


a.  to  seek  out  equipment  shortcomings  or  deficiencies  and 

to  determine  the  degree  of  additional  equipment  redundancy 
needed  for  future  testing; 

b.  to  gain  preliminary  insight  into  the  equipment  accuracy 
and  behaviour  in  the  overall  system;  and 

c.  to  ascertain  the  degree  of  control  necessary  to  hold 
the  control  parameters  within  test  code  bounds. 

3 . 0  REPORT  OBJECT 1 VE_S_ AND  SCOPE 

This  report  has  been  written  to  document  progress  of  intramural 
research  at  DREO  and  provide  the  necessary  background  for  DREO  and 
contract  personnel  working  in  the  Fuels/Powerplant  Technical  Subprogram. 

It  does  not  attempt  to  describe  equipment  choice  rationale  or  system 
design  parameters  but  describes  FEES,  the  tests  used  to  fulfill  the 
above  objectives,  the  results  obtained,  the  conclusions  and  recommendations 
pertinent  to  future  testing. 

4 . 0  THE  UTILIZATION  OF  FEES  AS  A  RESEARCH  TOOL 

FEES  has  been  designed  to  identify  small  differences  in  engine 
performance  and  operability  corresponding  to  either  subtle  or  major 
changes  in  fuel  chemistry.  Initially  only  a  comparison  of  the  gross 
parameters,  brake  power,  torque,  brake  specific  fuel  consumption,  volu¬ 
metric  efficiency,  air-fuel  ratio  and  brake  thermal  efficiency  will  be 
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The  most  >>asi<  part  01  tin*  c-stem  is  th-.  test  bed,  a  reinforced 
concrete  slab  on  top  of  a  db  mm  rubber  pad.  i  lu  let  bed  essentially 
rigidizes  tiie  dvn.im* I  ct— eng  i  iu*  i  al.ci  1  ofircol  ion  and  tin*  rub!>or  pad  and 
slab  dampens  v  ibr.H  ion  to  tin-  bui  ldinj:.  Hr-  contiol  room  is  used  to 
house  the  engine  --dvnar.tor.et  of  controls,  record  i  nr  instrumentation  and 
sensor  readouts.  Figures  da  and  db  sbow  ;  lie  eng i no -dynamometer  controller, 
strip  chart  recorder,  sensor  displ  readouts  and  a  data  acquisition 
system . 
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On  both  ends  of  the  toryuemelcr ,  flexible  rubber  couplings 
accommodate  slight  misalignment  and  prevent  shock  loads  being  transmitted 
from  rhe  engine  to  the  dynamometer.  In  order  to  reduce  power  dissipation 
due  to  misalignment,  an  alignment  jig  was  developed  to  permit  periodic 
checks  of  engine/torquemeter  alignment  with  the  dynamometer. 

Figure  4  a  schematic  drawing  of  the  entire  system,  itemizes 
measurement  locations  and  lists  recording  or  readout  devices.  As  can  be 
observed,  a  number  of  parameters  ate  measured  twice  to  improve  the 
experimental  reliability. 


Figures  5a  and  5b  indicate  the  sensor  locations  of  the  schematic 
relative  to  the  engine-dynamometer  position  on  the  test  bed.  In  greater 
detail,  Appendix  A  lists  each  of  the  measured  parameters,  the  symbol 
designated  to  it,  and  the  .on.-  i <: g  measurement  device. 
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To  permit  the  initiation  of  testing  as  quickly  as  possible,  FEES 
had  its  carefully  documented  circuitry  wired  in  a  temporary  fashion.  This 
also  facilitated  the  inevitable  changes  that  this  study  would  recommend. 

To  familiarize  the  reader  with  ti.e  KM914A  engine,  specification 
sheets  have  been  inserted  as  Appendix  B.  As  can  be  seen  from  the 
Appendix  B  performance  curves,  the  engine  is  capable  of  producing  12  kW 
at  4500  RPM.  Its  peak  torque  26  NM  occurs  at  3750  RPM  and  its  specific 
fuel  consumption  is  a  minimum  of  450  g/kWh  at  4100  RPM.  The  engine  is  a 
single  rotor,  carburetted  (Tilotson)  spark  ignited  rotary  combustion  prime 
mover.  Starting  is  accompl i.shed  by  using  the  D.C.  dynamometer  in  the 
motoring  mode. 

A  refrigeration  system  capable  of  temperatures  down  to  -40°C  has 
been  fabricated.  The  use  of  the  D.C.  dynamometer  with  the  engine  in  a 
refrigerated  enclosure  provides  suitable  cold  cranking  capability.  The 
proposed  addition  of  combustion  air  and  cooling  air  temperature-humidity 
control  to  -40°C  should  also  permit  low  temperature  engine  operation  for 
fuels  testing  of  research  type  single  cylinder  engines. 

6 . 0  TEST  PROGRAMME 

It  was  decided  to  conduct  tests  at  three  different  throttle 
openings  50%,  75%  and  100%.  This  was  accomplished  in  a  reproducible 
manner  using  a  throttle  actuator  and  electronic  control  device.  At  each 
throttle  opening  the  dynamometer  load  was  set  to  determine  the  operating 
RPM.  In  this  fashion  the  RPM  was  varied  between  2000  and  4000  RPM.  For 
a  particular  test  run  at  a  fixed  throttle  opening  and  adjusted  load 
condition,  the  engine  was  operated  for  approximately  5  to  15  minutes  in 
a  pretest  period  until  temperature  stabilization  occurred  at  the  exhaust 
manifold  thermocouple.  Once  temperature  stability  was  present  the  test 
run  was  begun  and  readings  were  recorded  for  a  know  time  period  of 
approximately  four  minutes.  After  this  time  period  elapsed,  exhaust 
manifold  temperature  was  again  observed.  If  this  continued  to  be  stable, 
results  were  calculated  immediately  and  torque  and  brake  power  plotted 
to  compare  with  previous  test  runs.  During  this  post  test  period  of 
two  to  three  minutes  other  results  were  fed  into  the  DREO  site  main  frame 
computer.  If  the  plotted  data  compared  favourably  and  the  temperature  of 
the  exhaust  gases  and  other  temperatures  remained  stable,  a  new  load 
setting  was  effected  and  pretest  checking  for  a  new  test  run  was 
commenced. 

To  minimize  errors  and  omissions,  with  the  large  number  of  sensors 
and  variety  of  recording  or  readout  devices,  elaborate  check  lists  were 
devised.  Check  lists  at  the  beginning  and  end  of  the  testing  day  ensured 
proper  calibration  of  all  equipment.  Pretest  check  lists  ensured  all 
previous  data  was  recorded  correctly  and  that  the  next  data  would  be 
taken  in  the  correct  format.  Due  to  the  thoroughness  of  check  lists,  no 
errors  or  omissions  resulted  during  the  14  tests  for  this  study. 
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The  results  of  a  series  of  tests  for  a  particular  throttle  opening 
led  to  the  calculation  of  the  performance  parameters.  These  performance 
parameters  outlined  in  the  following  section  are  treated  in  greater  detail 
and  their  interrelationships  discussed  in  Taylor  [1],  Obert  [2]  and 
Lichty  [3].  From  the  performance  parameters,  performance  curves  can  be 
plotted  which  form  the  basis  for  comparisons  between  test  fuels  and  a 
baseline  or  reference  fuel.  In  addition,  the  engine  can  be  observed  for 
its  operability  characteristics ,  i.e.  how  effectively  or  smoothly  it 
operates  with  a  particular  test  fuel.  (This  observation  can  be  quantified 
by  utilizing  accelerometers  mounted  in  various  engine  locations  to  measure 
combustion  roughness  or  vibrational  signature.) 


7 . 0  GROSS  PERFORMANCE  PARAMETER  EQUATIONS 

The  gross  performance  parameters  are  calculated  from  the  measured 
parameters  described  in  Table  A1  of  Appendix  A  in  the  following  manner. 
Engine  steady  state  output  torque  (T)  is  read  directly  from  the  torque- 
meter  in  (in/lbs)  and  converted  to  (NM)  or  calculated  by 


T  =  F  x  r  (NM)  (1) 

s 

where  Fs  is  the  dynamometer  reaction  force  and  r  is  the  radius  measured 
from  the  dynamometer  centre  line  to  the  point  of  force  measurement.  The 
brake  power  (bp)  is  calculated  using 

bp  =  (Watts)  (2) 

where  N  =  Engine  Speed  (RPM) 

T  =  Torque  (NM) . 


The  brake  specific  fuel  consumption  (bsfc)  is  the  mass  flow  rate 
of  the  fuel  (mj )  divided  by  the  brake  power  (bp) . 
m. 


BSFC  = 


bp 


(3) 


Using  the  lower  calorific  value  of  the  fuel  (LCV)  determined  from 
API  gravity  and  aniline  point  for  the  ASTM  (D1405-64)  test  code,  the 
brake  thermal  efficiency  (ng^)  can  be  calculated  from 


bp’ _ 

n6T  m  x  LCV 


w 


Another  performance  parameter  is  the  cylinder  mean  effective 
pressure  which  can  be  expressed  as  the  brake  mean  effective  pressure 
(bmep) . 


The  torquemeter  used  also  had  the  ability  to  measure  instantaneous  torque 
of  the  engine  output  shaft  for  monitoring  combustion  space  energy  transfer. 
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bmep  =  -  -  b^-  - 
DT  x  N 


(5) 


where  (D.p)  is  the  total  displacement  and  (N)  is  the  engine  rotational 
speed  in  RPM.  The  bmep  may  be  thought  of  as  that  mean  effective  pressure 
acting  on  the  pistons  which  would  give  the  measured  bp  if  the  engine  were 
frictionless. 

For  air  breathing  engines  air  capacity  (ma)  is  an  important  para¬ 
meter.  Its  effect  on  indicated  power  (IP)  can  be  seen  in  the  following 
equation. 


IP 


J"a  <FVth> 


'6) 


where:  IP  -  power  developed; 

J  =  mechanical  equivalent  of  heat; 

ma  =  mass  flow  rate  of  dry  air; 

Qc  =  heat  of  combustion; 
rijH  =  indicated  thermal  efficiency;  and 
F  =  fuel/air  ratio. 

For  SI  engines,  the  indicated  power  remains  proportional  to  air  capacity, 
provided  there  is  no  change  in  fuel/air  ratio  or  compression  ratio  and  no 
departure  from  optimum  spark  timing.  For  these  conditions,  indicated 
thermal  efficiency  remains  reasonably  constant  and  this  proportionality 
between  IP  and  ma  has  been  found  to  hold  for  many  types  of  SI  engines. 

To  further  examine  the  aspiration  capability  of  an  engine,  it  is 
convenient  to  use  an  expression  independent  of  cylinder  size.  Such  an 
expression  is  the  volumetric  efficiency.  For  ideal  engines  volumetric 
efficiency  (hv)  may  be  defined  as 


%  = 


(Vi-V7)Ptr 


(7) 


where  m^  is  the  mass  of  fresh  mixture  of  moist  air  (m)  supplied  and  is 
the  density  at  the  pressure  pm  and  temperature  Tn].  The  quantity 
(V^-V2)pm  is  the  mass  of  fresh  mixture  which  would  just  fill  the  piston 
displacement  volume  at  the  density  of  the  inlet  system.  The  volume 
is  at  the  piston  bottom  dead  center  (bdc)  (position  -x^)  and  the  volume  V 
is  at  top  dead  center  (tdc)  (position  -X2) . 


For  real  engines  the  volumetric  efficiency  can  be  expressed  as 


2m 
_  m 

nv  -  NVdJPn 


(8) 


where  =  mass  flow  rate  of  fresh  mixture  (m) ; 

N  =  RPM; 

Vd  =  total  displacement  volume  of  engine; 

Vd  =  (x^-X2)Ap  where  Ap  =  bore  area  of  piston;  and 
pm  =  inlet  density  of  fresh  mixture  (m) . 
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The  factor  2  in  equation  (8)  arises  from  the  engine  combustion  cycle.  The 
2  identifies  a  four  stroke  cycle  where  there  is  one  complete  cycle  every 
two  crankshaft  revolutions.  The  equation  would  be  without  the  2  for  a  two 
stroke  combustion  cycle. 

For  rotary  combustion  engines  (RCE)  there  has  been  extensive 
disagreement  as  to  what  constitutes  the  effective  displacement  volume  and 
whether  the  engine  operates  in  a  two  or  four  stroke  cycle.  Norbye  [4] 
indicates  that  RCE ' s  are  now  classified  as  following  the  four  stroke 
cycle  and  that  the  Commission  Sportive  Internationale  (CSI)  has  adopted  a 
formula  which  rates  RCE's  displacement  at  twice  the  combustion  chamber 
volume  multiplied  by  the  number  of  rotors.  This  eliminates  the  2  in 
equation  (8)  and  implies  RCE's  breathe  much  like  a  two  stroke  cycle  power- 
plant. 


For  equation  (8),  it  is  necessary  to  define  inlet  density  as  the 
density  of  fresh  mixture  in  or  near  the  inlet  port.  The  resulting  volu¬ 
metric  efficiency  measures  the  pumping  performance  of  cylinder  and  valves 
alone.  Since  it  is  not  always  possible  or  convenient  to  measure  pm  at 
the  inlet  port,  density  is  usually  measured  at  the  engine  air  intake. 

The  volumetric  efficiency  based  on  this  measurement  location  is  called 
the  overall  volumetric  efficiency.  For  unsupercharged  engines  with  small 
pressure  and  temperature  changes  in  air  cleaner,  carburettor  and  inlet 
manifold  geometry  the  overall  volumetric  efficiency  does  not  differ 
greatly  from  the  port  volumetric  efficiency. 


The  fuel/air  ratio  (F)  is  calculated  from 


where  iiif  =  mass  fuel  flow  rate 

n>a  =  PaQa 

where  Qa  =  volumetric  flow  rate  of  dry  air 
pa  =  density  of  dry  air. 


(9) 


The  air/fuel  ratio  (A/F)  is  simply 


(10) 


It  is  appropriate  to  introduce  the  term  equivalence  ratio  (<P) ,  which  is  the 
quantities  of  equation  (9)  and  (10)  normalized  by  the  stoichiometric  or 
chemically  correct  fuel/air  ratio 

.  _  F(ac tual)  _  A/F( s toich)  Mi} 

F(stoich)  A/F(actual) 


For  octane  (CgH^g)  F(stoich)  is  usually  given  as  0.0668. 

These  performance  equations  have  been  used  in  a  sample  calculation 
in  Appendix  D  for  test  run  No.  RC30.  The  test  conditions  for  test  run 
No.  RC30  are  tabulated  in  Appendix  C,  Tables  Cl  and  C2  as  are  the  data 
of  the  earlier  tests  RC16  to  RC29. 
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8.0  PERFORMANCE  RESULTS 


TABLE  2 

Test  Data  For  API  Gravity  And  Calculated  Higher  Heating  Value 


RUN 

NO. 

Z 

(%) 

API 

'  higher' 

HEATING 

VALUE 

(BTU/LB) 

RC  (16-20) 

50 

65.8 

20552 

RC  (21-25) 

75 

64.8 

20512 

RC  (26-30) 

100 

62.0 

20400 

Test  Altitude  = 

232  feet 

above  sea 

level 

The  performance  equation  data  has  been  tabulated  in  Tables  1  and  2 
and  plotted  against  engine  speed  in  the  graphs  of  Figures  6  to  12.  From 
the  torque  values  of  Figures  6a,  a  smooth  curve  is  present  for  the  1 00% 
throttle  opening  case. 


E 

z 


ENGINE  SPEEO  (N(  rpm 


Figure  6a  -  Torque  vs  Engine  Speed  for  the  Sachs 
KM914A  Engine 
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13. 


For  the  75%  and  50%  part  throttle  opening  curves,  scatter  in  data  was 
primarily  due  to  difficulties  in  reading  the  dynamometer  spring  balance. 
Severe  vibrations  at  certain  frequencies  caused  erratic  needle  deflection, 
that  was  difficult  to  read  accurately.  These  vibrations  were  observed 
to  be  a  function  of  engine  speed  and  load.  Fortunately  excessive 
vibration  did  not  occur  for  any  of  the  runs  of  primary  importance  at  100% 
throttle  opening.  It  can  be  seen  that  smooth  curves  exist  for  all  the 
air  mass  flow  rate  and  fuel  mass  flow  rate  data  of  Figures  6b  and  6c 
respectively.  This  indicates  that  the  scatter  present  in  Figure  6a  is 
not  engine  related  but  dynamometer  related. 


ENGINE  SPEED  <Nl  *p>n 


Figure  Cb  - 


Intake  Mr  Haas 
Engine  Speed 


Flow  Rate 


On 


it 
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Figure  6a  -  Fuel  Hass  Flow  Fate  vs  Engine  Freed 

The  electric  dynamometer  proved  to  be  difficult  to  control 
throughout  the  tests.  The  D.C.  power  to  the  dynamometer  field  coils  was 
controlled  by  a  potentiometer  which  was  continually  adjusted  by  hand  to 
keep  a  standing  pressure  wave  form  generated  by  the  combustion  space 
pressure  transducer  centred  on  the  oscilloscope  screen  for  a  given 
engine  rotational  frequency.  This  yielded  speed  control  to  within  15  RPM 
and  load  control  for  dynamometer  torque  of  4  NM.  Engine  torque  is  usually 
measured  by  the  shaft  torque  meter  shown  in  Figure  3,  but  due  to  brush 
problems  with  the  torque  meter,  the  backup  spring  balance  was  used  for 
all  of  these  tests. 

From  equation  (2)  the  brake  power  (bp)  values  of  Table  1  were 
generated  and  the  smooth  curves  of  Figure  7  were  plotted.  Despite  minor 
scatter  at  low  speeds,  the  crude  speed/load  control  and  the  inaccurate 
spring  balance,  the  bp  results  were  encouraging.  These  results  suggested 
that  bp  accuracy  and  repeatability  within  1.5%  for  fuels  testing  would  be 
possible  with  improved  dynamometer  control  and  use  of  the  torque  meter. 
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ENGINE  SPEED  INI  rpm 


Tnem.x'l  Efficiency  vo  Engine  Speed 

I  row  the  brake  thermal  efficiency  curves  of  Figure  9,  it  is  evident 
that  this  particular  rotary  combustion  engine  is  highly  inefficient.  From 
equation  (4)  and  Figures  6c  and  7,  the  high  fuel  consumption  at  3500  RPM 
yielded  a  pronounced  drop  in  efficiency  in  Figure  9.  It  is  suspected  that 
the  poor  performance  of  the  engine  in  terms  of  BSFC  and  n,:-f  was  due  to 
compression-combustion-expansion  leakage  past  the  rotor  face  and  apex 
seals.  This  leakage  should  result  in  a  lower  brake  power  which  was  not 
the  case.  The  brake  power  in  Figure  7  compared  favourably  with  the 
manufacturer's  ratings.  The  reason  for  this  anomaly  is  related  to  exhaust 
back  pressure.  Had  the  engine  been  tested  with  higher  back  pressure 
equivalent  to  the  manufacturer’s  muffler,  the  output  would  have  been 
considerably  lower  as  would  be  expected  with  leaking  seals. 
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ENGINE  SPEED  N  ipn- 

Figure  10  -  Brake  Mean  Effective  Pressure  vs  Engine  Speed 

The  brake  mean  effective  pressure  of  .55  mPa  in  Figure  10  for  Z 
100%  at  4000  RPM  is  about  2%  lower  than  a  performance  map  value  in 
reference  [4].  As  expected,  the  general  shape  of  the  BMEP  curves  is  the 
same  as  the  brake  power  curves  in  Figure  7  since  brake  power  is  used  to 
calculate  bmep  of  equation  (5). 


Figure  11  -  Volumetric  Efficiency  vs  Engine  Speed 
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From  Figure  11  it  can  be  seen  that,  at  full  throttle  opening,  the 
volumetric  efficiency  increases  linearly  with  engine  speed  and  reaches  a 
value  in  excess  of  100%.  This  is  typical  of  rotary  combustion  engines 
as  they  breathe  more  efficiently  than  do  reciprocating  engines.  One 
might  think  of  a  rotary  engine  as  a  type  of  engine  that  lies  somewhere 
between  a  reciprocating  and  gas  turbine  powerplant. 

The  part  throttle  curves  of  Figure  11  can  best  be  explained  by 
the  air/fuel  ratio  curves  of  Figure  12.  The  full  load  curve  is  smooth 
without  discontinuities,  whereas  the  75%,  and  50%  curves  show  the  effects 
of  throttling  the  air  flow  before  the  combustion  process.  This  process 
usually  leads  to  so  called  part  throttle  performance  flat  spots  that  are 
evident  at  2500  and  3500  RPM. 


ENGINE  SPEED  iN  -pm 


Fi-r.iVc  -  A:i‘/F:u'!  Fatee  Fi  cc.i 

The  overall  performance  data,  as  tabulated  in  Tables  1,  2,  Cl  and 
C2,  are  in  good  agreement  with  the  values  plotted  in  the  curves  of 
Figures  6  to  12.  The  notable  exception  is  the  spring  balance  reading. 

It  is  planned  to  replace  this  device  with  torque  strain  gauge  transducers 
in  future  testing.  It  was  also  determined  through  these  tests  that  the 
sampling  rate  of  the  data  acquisition  system  used  for  temperature 
measurement  could  not  be  increased  to  accommodate  transient  pressure 
time  information. 

9 . 0  CONCLUSIONS  AND  RECOMMENDATIONS 

The  major  equipment  short-coming  uncovered  as  a  result  of  these 
tests  is  related  to  the  measurement  of  torque.  It  is  recommended  that  the 
dynamometer  spring  balance  reading  be  replaced  by  a  strain  guage  brushless 
torque  transducer.  Such  a  transducer  must  be  capable  of  reading  both 
steady  state  and  dynamic  torque.  The  second  equipment  short-coming 
concerns  the  data  acquisition  system.  Although  suitable  for  temperature 
measurement,  the  system  was  shown  to  have  a  sampling  rate  inadequate  for 
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monitoring  pressure-time  or  other  engine  transients.  It  is  recommended 
that  a  Tektronix  Digital  Analyzer  system  be  considered  for  monitoring 
transient  engine  phenomena. 

All  other  equipment  performed  as  expected  in  either  a  primary  or 
backup  measurement  role.  It  is  believed  that  performance  reading  for 
primary  measurements  will  be  achievable  to  within  1  1/2  %.  This  will 
enable  the  system  to  detect  small  changes  in  performance  as  a  result  of 
subtle  changes  in  fuel  chemistry.  It  is  recommended  that  all  equipment 
now  be  hardwired  into  the  system  and  coupled  to  a  Tektronix  Digital 
Analyzer  and  Fluke  data  logger. 

In  order  to  control  the  engine  load  provided  by  the  D.C.  dynamometer 
during  testing  a  different  approach  than  the  manually  controlled  field  coil 
resistor  will  have  to  be  adopted.  Control  to  within  +  10  RPM  is  necessary 
to  comply  with  SAE  specifications.  It  is  recommended  that  an  electronic 
controller  be  built  or  purchased.  Such  a  controller  will  control  the 
field  current  in  response  to  a  Hall  effect  speed  sensor  with  a  60  tooth 
gear.  This  should  provide  adequate  control  for  upcoming  diesel  engine 
fuels  testing. 
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APPENDIX  A 


PARAMETER  DESCJR I PT_I ON,  SYMBOL  AND 
CORRE ;s PONDING  MEASUREMENT  DEVICE 
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TABLE  A1 


Parameter  Description  and  Measurement  Device  List 


PARAMETER 

NAME 

UNITS 

PARAMETER  DESCRIPTION 

MEASUREMENT  DEVICE 

Pb 

In  Hg 

Atmospheric  pressure 

Barometer 

td 

°F 

Temperature  Dry  Bulb 

Thermometer 

T 

W 

°F 

Temperature  Web  Bulb 

Sling  Psychrometer 

0 

% 

Relative  Humidity 

Psychrometric  chart 

mm  Hg 

Vapor  Pressure 

Vapor  Pressure  chart 

- 

Correction  factor  for 
Humidity 

Humidity  Correction 
chart 

API 

- 

Fuel  Specific  Gravity 

Hydrometer 

°NOZ 

in 

Diameter  of  Air  Drum 

Nozzle 

Pre-calibrated 

Z 

o 

Throttle  Opening 

Throttle  Activator 

L 

- 

Load 

Dynamometer  Load 
Position  Knob 

I 

Amps 

Load  Current 

Dynamometer  Guage 

V 

Volts 

Load  Voltage 

Dynamometer  Guage 

T 

in  lbs 

Engine  Torque 

Torquemeter  (Lebow) 

Qf 

cc  sec~l 

Fuel  Flow  rate 

Positive  Displ  Flow 
Meter  (Fluidvne) 

QA 

ft^  tnin~l 

Air  Flow  rate 

Vane  Type  Flowmeter  j 

(AutoTronics) 

T 

min&secs 

Elasped  time  for  run 

Digital  Clock 

VT 

cc 

Total  Volume  of  Fluid 

Positive  Displ  Flow¬ 
meter  (Fluidyne) 

0 

c 

Radians 

Crank  Angle  Position 

Hall  Effect  Velocity 
Sensor  (Airpax) 

0 

Radians 

Power  Take  off  Shaft 

Hall  Effect  Velocity 

s 

(PTD)  Position 

Sensor  (Airpax) 

nca 

rpm 

Revolutions  per  minute 
of  Crankshaft 

Hall  Effect  Velocity 
Sensor  (Airpax) 

P 

C 

psi 

Cylinder  Pressure 

Piezo  Electric  Pressure 
Transducer  (Kistler) 

PA 

in  H20 

Intake  Air  Pressure 
across  orifice 

Intake  Air  Flow  Measur¬ 
ing  (Go  Power) 
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APPENDIX  B 


MANUFACTURER  SPECIFICATIONS  FOR  THE  KM914A  SACHS, 
SPARK  IGNITED,  ROTARY  COMBUSTION  ENGINE 
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APPENDIX  C 


COMPUTER  PRINTOUT  FOR  TEST  RUN  #RC3Q 
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TABLE  Cl 

Computer  Printout  for  Test  Run  #RC30 


RUN  BENGINE1 
PI  ASSOCIATED 

GROUP  17  OF  : LIB.: SYS  ASSOCIATED  FOR  9INPUTU 
*  *  ALLOCATION  SUMMARY  *  * 

PROTECTION  LOCATION  PAGES 

DATA (00)  A000  3 

PROCEDURE  (01)  A800  7 

DCB  (10)  A600  1 

FILE  NAME  =  RC30 


ATM  Pressure 

= 

30.120 

IN  HG 

DPB 

Temp  Dry  Bulb 

= 

53.000 

DEC.-F 

DTDF 

Temp  Wet  Bulb 

= 

40.000 

DEG-F 

DTWF 

Vapor  Pressure 

= 

.000 

MM  HG 

DE 

Relative  Humidity 

= 

25.000 

% 

DPHI 

Corr.  Factor  For  Humidity 

= 

1.000 

DHC 

API  Number 

= 

62.000 

DEG 

DAP  I 

Diameter  of  Air  Drum 

= 

.750 

IN 

DNOZ 

Throttle  Opening 

= 

300.000 

/300 

DZ 

Load 

= 

15.000 

DL 

Current 

= 

54.000 

AMP 

DI 

Voltage 

= 

177.000 

VOLTS 

DV 

Torque  (Lebow) 

= 

.000 

IN-LBF 

DTLS 

Fuel  Flow 

= 

2.260 

CC/SEC 

GFUEL 

Air  Flow 

= 

796.000 

HZ 

DAA 

Time 

= 

400.000 

MIN-SEC 

DFT 

Volume  of  Fuel 

= 

546.420 

CC 

DFVF 

Period  of  Crank 

= 

.000 

SEC 

TETAC 

Period  of  PTO  Shaft 

= 

.000 

SEC 

TETAS 

RPM  (Crank) ,  Analog 

= 

4000.000 

RPM 

NCA 

Cylinder  Pressure 

.000 

PSI 

PCYL 

Air  Drum  Pressure 

= 

.000 

IN  H2O 

DAM 

RPM  (Crank) ,  Photo 

= 

4000.000 

RPM 

DN1 

Force  (Dynamo) 

= 

55.000 

NEWTONS 

DSBR 

RPM  (PTO),  Dynamo 

= 

2274.000 

RPM 

DN2 

Mass  of  Fuel 

= 

.393 

KG 

DFMB 

No.  of  Revolutions 

= 

4528.000 

REV 

DREV 

TEMPERATURES: 

Initial 

(1) 

679.800 

(2) 

659.100 

(3) 

13.300 

(A) 

10.300 

(5) 

8.100  C 

Final 

683.2 

661.1 

14.2 

10.1 

8.7  C 
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TABLE  C2 

Computer  Printout  for  Test  Run  //RC30 

FILE  NAME  =  RC  30 


Relative  Humidity  =  2.625/10.289  =  25.5% 


Barometric  Pressure 

30.12 

in  of  Hg 

765.05 

mm  of  Hg 

Dry  Bulb  Temperature 

11.67 

C 

53.00 

F 

Wet  Bulb  Temperature 

4.44 

C 

40.00 

F 

Relative  Humidity 

25.5  % 

Altitude  Above  SDA  Level 

232.00 

feet 

70.71 

M 

Vapor  Pressure 

2.63 

mm 

70.71 

m 

Load 

15 

Current  Dynamometer 

54.00 

ADC 

Voltage  Dynamometer 

177.00 

UDC 

Power  Dynamometer 

12.82 

hp 

9.56 

kw 

Density  of  Fuel 

.731  C/CC 

Density  of  Air 

1.247  G/L 

Torque  (Lebow  Sensor) 

.00 

in-lbs 

.00 

N-M 

Torque  (Spring  Balance) 

414.21 

in-lbs 

46.80 

N-M 

Fuel  Flow  Rate  (Fluidyne) 

Volume : 

2.28 

CC/sec 

Mass : 

1.66 

G/ sec 

.29 

CF/hr 

13.21 

lbs/hr 

Fuel  Flow  Rate  (Balance) 

Volume : 

2.24 

CC/sec 

Mass : 

1.64 

G/sec 

.28 

CF/hr 

13.00 

lbs/hr 

Air  Flow  Rate  (Autotronics) 

Volume: 

21.00 

L/sec 

Mass : 

26.18 

G/sec 

2670.25 

CF/hr 

207.79 

lbs/hr 

Air  Flow  Rate  (Manometer) 

Volume : 

.95 

L/sec 

Mass : 

1.18 

G/ sec 

120.42 

CF/hr 

9.37 

lbs/hr 

Air/Fuel  Ratio 

15.73 

15.99 

.71 

.72 

Brake  Power 

14.95 

hp 

11.15 

kw 

.00 

hp 

.00 

kw 

Brake  Specific  Fuel  Consumption 

.87 

lbs/hp-hr 

.53 

kg/kw-hr 

#### 

lbs/hp-hr 

IHHHt 

kg/kw-hr 

.88 

lbs/hp-hr 

.54 

kg/kw-hr 

////## 

lbs/hp-hr 

#### 

kg/kw-hr 

Brake  Mean  Effective  Pressure 

80.45 

lbs/sq  in 

.00 

lbs/sq  in 

Thermal  Efficiency 

14.35 

% 

.00 

% 

14.12 

% 

.00 

% 

Volumetric  Efficiency 

104.49 

% 

4.71 

X 
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Table  C2  (Cont'd) 


Power  Distribution:  IP  =  BP+FP 
20.196  =  .950  +  19.179  hp 

Corrected  We  Obtain: 

Power  Distribution:  IP  =  BP+FP 
20.196  =  14.952  +  5.244  hp 

Corrected  We  Obtain: 

19.2  =  13.9  +  5.2 

TC  =  386.0193  in-LBF 

BMEPC  =  74.9725  psi 
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APPENDIX  D 

GROSS  PERFORMANCE  PARAMETER  CALCULATIONS 
FOR  TEST  RUN  //PC  30 
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The  sample  calculations  performed  for  this  Appendix  are  for  test 
run  #RC30.  The  test  conditions  are  listed  as  the  input  to  a  computer 
program  in  Appendix  C,  Table  Cl.  The  results  of  the  program  are  listed 
in  Table  C2.  These  results  can  be  compared  with  the  sample  calculations 
for  this  Appendix.  For  brevity,  only  one  set  of  calculations  has  been 
submitted.  These  calculations  use  data  from  the  instrumentation  considered 
to  be  the  most  accurate.  For  instance,  calculations  of  the  air  flow 
are  based  on  the  air  turbine  rather  than  the  less  accurate  pressure  drop 
across  the  orifice. 

The  gross  performance  equations  used  in  the  calculations  are 
discussed  in  Section  7.0  in  the  main  body  of  the  report. 


Engine 

Shaft  Rotation 


Figure  D1  -  Dynamometer  Torque  Measurement 

The  torque  in  equation  (1) 

T  =  F  x  r  (1) 

s 

is  calculated  by  finding  the  true  rotational  force  opposing  engine  shaft 
rotation.  This  can  be  found  by  merely  subtracting  the  dynamometer  spring 
indicator  dial  reading  (Fp)  from  the  equilibrium  condition  of  zero  torque 
i.e.  250N  of  the  counter  weight. 

F  =  250N  -  F 
s  D 

From  the  computer  printout  in  Appendix  C  for  test  RC30,  the  input  section. 
Table  Cl  lists  the  dynamometer  dial  indicator  force  reading  (Fp)  as  55N. 

Fg  =  (250  -  55)N  =  195N 

The  force  arm  radius  (r)  in  equation  (1)  and  Figure  D1  is  0.240m. 
Thus  torque  is 
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T  =  195N  x  0.24m  39.37  in  =  1  in 

=  46.8  Nm  (414.2  in  lb)  @  2275  RPM  4.4482  N  =  1  lbf 

The  brake  power  (bp)  can  be  calculated  by  substituting  this  value  of  T 
and  the  Dynamometer  RPM  from  Table  Cl  in  equation  (2).  The  Dynamometer 
RPM  of  2275  is  used  rather  than  the  engine  RPM  of  4000  because  the  torque 
is  measured  on  the  dynamometer  side  of  the  1.758/1  gear  reduction  unit. 
The  torque  before  gear  reduction  would  be 

T  =  46 . 8Nm/ 1.758  =  26.6  Nm  @  4000  RPM 


bp 


2ttNT 

60 


(2) 


2frrad  2275  rev  1  min  .  .  0  „  ,1  watt 

=  - - —  x  x  v ■-  —  ■  x  46.8  Nm  (-  -  -  - 

rev  min  60  sec  INm/sec 

2tt 

=  (■>—  x  2275  x  46.8)  watts 
60 

=  11,149.5  watts  or  11.15  Kw  (14,95  Hp) 


The  mass  fuel  flow  rate  (nif)  can  be  calculated  by  using  from  Table  Cl  the 
mass  of  fuel  393  Kg  over  the  test  period  of  4.0  min  (240  secs). 


393 

mf  "  240  sec 


1.64  g/sec 


The  brake  specific  fuel  consumption  (bsfc)  of  equation  (3)  can  now  be 
calculated 

,  - 

bsrc  =  b-p 

1.64  g.sec  3600  sec 
11.15  Kw  X  1  hour 


=  529.5  g/Kwh. 


It  is  n>  generally  accepted  that  the  lower  calorific  value  (LCV) 
be  used  in  determining  brake  thermal  efficiency  (equation  (4)).  However, 
for  the  purpose  of  these  shakedown  tests  it  was  considered  unnecessary  to 
carry  out  aniline  point  tests  for  each  fuel  batch  according  to  ASTM 
(D611-77).  Instead  of  using  the  conventional  means  of  calculating  LCV, 
ASTM  (D1405-69) ,  the  modified  Sherman -Kropf f  equations  were  used  to 
determine  the  higher  heating  value  (HHV). 

HHV  =  18,320  +  40  (API-10)*  in  —  (for  gasoline) 


*Journal  of  the  American  Chemical  Society,  Vol  30,  No  10,  p  1626. 
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This  is  an  approximate  HHV  equation  where  only  the  API  gravity  is  needed, 
ASTM  (D1298-67).  HHV  is  usually  about  6  to  9%  higher  than  the  LCV.  This 
will  tend  to  reduce  the  brake  thermal  efficiency.  As  Table  Cl  indicates 
the  API  gravity  was  experimentally  determined  in  accordance  with  ASTM 
(1298-67)  and  found  to  be  62. 

HHV  =  18,320  +  40  (62-10) 

=  18,320  +  2,080 

=  20,400  B-™  (47,405  |~) 

The  brake  thermal  efficiency  Ogj  tna>'  now  be  calculated  using  equation  (4) 
with  the  HHV  revision. 

bp 

nBT  mf  x  HHV 


_ 11.15  Kw _ 

1.64  g .sec  x  47,405  J/g 

_ li.15  x  l_o£ w _ 

1.64  g.sec  x  47,405  - 

S 

=  14.34% 

To  calculate  the  brake  mean  effective  pressure  (bmep)  from 
equation  (5)  an  equivalent  displacement  term  must  be  used.  For  the 
purposes  of  these  tests  the  manufacturers  equivalent  displacement  of 
303  cc  (18.49  in^)  per  chamber  will  be  used. 


bmep  = 


Total  Displacement  revolutions 
revolution  min 

11.15  x  103  N*m'sec  ^  _ _ 

. y  "  . . /  .  y 

303  cm  4000  rev  1  min  1  m 

rev  mm  60  sec  .  _6  3 

10  cm 


11.15  x  10 
303  x  A 


11.15  x  106 
’  20.4 

=  . 552  x  106  N 


=  .552  MPa  (79.33  psi) 
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In  order  to  calculate  the  overall  volumetric  efficiency  (r)y)  in 
accordance  with  equation  (8) ,  the  density  of  the  intake  air  mixture  of 
water  vapour  with  dry  air  must  ’  e  found.  As  mentioned  earlier  if  the 
volume  flow  rate  (Vm)  is  known  tor  equation  (8)  the  calculation  of  pm  is 
no  longer  necessary.  The  flow  device  used  for  this  study  had  the  option 
of  using  either  Vm  or  for  a  test  run  and  both  were  alternately  used. 
Since  the  reading  for  Mm  has  to  be  verified  pm  will  be  calculated. 

To  calculate  air-vapour  mixture  density  a  number  of  empirical 
equations  exist  and  for  atmospheric  pressures  between  .9  and  1.1  bar  a 
psychrometric  chart  is  accurate  to  ±  3%.  For  completeness  and  benchmark 
purposes  a  fundamental  approach  has  been  used  here.  English  units  have 
been  used  for  the  calculations  with  final  air  densities  being  converted 
to  metric  units.  For  a  more  complete  discussion  of  the  following  subject 
it  is  suggested  the  reader  refer  to  any  classical  thermodynamics  text. 

A  temperature  entropy  (t-s)  diagram  can  be  drawn  as  in  Figure  D1 
where  the  environmental  conditions  for  test  RC30  are  shown.  The  dry  bulb 
temperature  (t^  =  53°F)  and  the  wet  bulb  temperature  (tw  =  40°F)  are 
located  for  a  barometric  pressure  (pni  =  30.12  in  of  Hg) . 


>-  I 
3  ! 

•u 

U 

cy 

a. 

E 

<U 


td=53°F  b 


A 


t  =40°F 
w 


Dew  point  temperature  t 


dp 


Entropy  -  s 

Figure  D2  -  Temperature  Entropy  Diagram  for  Test  RC30 

For  an  adiabatic  saturation,  irreversible  process  an  energy 
balance  for  Figure  D3  yields 


h  ,  +  w ,  h  ,  +  (to  -  u),)hc  =  h  +  w  h 
ad  d  vd  w  d  fw  aw  w  vw 


(1-D) 


,  lb  of  vapour  ,  . 

where  u>,  =  qr — c  f  —  -  and  evaporation  =  to  -to, 

d  lb  of  dry  air  w  d 
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43. 


Figure  D2  -  Adiabatic  Saturation  Process 


The  enthalpy  of  evaporation  symbol 

h  ,  -  h  -  hP 

fgw  vw  fw 

for  the  exit  state  equation  1-D  may  be  solved  for  the  humidity  ratio  (tod) 
of  the  entering  air. 


0) 


d 


^w^fgw^aw  ^ad  ,1b  of  vapour  . 

h  ,-h,  'lb  of  dry  air 

vd  fw 


(2-D) 


where 

h  -h  .  * 

C  At 

aw  ad 

P  a 

and 

t  *  t  - 

t.,  when  C  ®  .24 

a  w 

d  p 

From  thermodynamic  steam  tables  for 
in  Hg 


BTU/lb-°R 

tw  =  40°F,  td  =  33°F  and  pm  =  30.12 


h, 

fgw 


1071 


BTU 

lb 


1084.7 


BTU 

lb 


fw 


8.027 


BTU 

lb 


At  t  =  40°F  the  partial  pressure  of  the  steam  is  pvw  =  .12163  psia  and 
fromWsteam  tables  Vg  =  2445.8  ft3/lb  or  pw  *  1/2445.8  =  4.089  x  10-^ 
lb/ft3 


Pm 


30.12  in  of  Hg  x 


.491  psia 
1  in  of  Hg 


14.78  psia 
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n  1  !  ',  !  ,i-.ia 


•  1 1-  ;i  I  .*  t  (,-J  us  i  a,  tin:  perfect  gas  assumption 


t  '  i  •  )  tii  li:  ,  (  in'  i l  " ) 

s  s.  3  (  V.".'  K;  it  lli-  <  lb°K>  1  °R 


-  .  0 .  v  2  5  1  b  /  r  t 

■  i  t  v  ratio  after  sat  :  .  .  u.  into  ■■  now  be  found . 


vv:  _  4 .  lifer1  x  i  i 
j  .07  925 


j  .  1  59  10 


3  ! 1 


Ti.e  ;ir.  ••  h.}i  tel  at  i. 'ns.-,  ip  tor  equation  (2-0)  may  be  found  where  the  state 
1  .subscript  replaces  subscript  J. 

huv  -  haj  -Cp(t  -  -0.14(53-40) 

The  iiuniuitv  ratio  u-j  of  the  original  air  from  equation  ( 2— D )  can  now  be 
calc.ul  at. cd 


_  Ww  "fgw  +  '' aw  !lai 

U1  ‘  '  hv'l  -  hfw  ’ 

(5.159  x  10~3)  (1071)  -  (0.24_(53-40i_) 
1084.7  -  s’.  02  7 


=  2.2  34  x  10' 


.  •}  lb  v 
lb  d,. 


The  partial  pressure  of  Liu.-  vapour  is  from  Carrier's  equation 

1' 


(3-D) 


1 


-  ls 


Pvp  - 


(2.  2  34  x  10  )  (  I  ..  78 ) 

(0.022)  +  (2.234  x  10~3) 


5.289  x  10  “  psia. 


At  t^  =  53°F  the  f 


saturation  pressure  from  steam  tables  is  pv^  =  .19883 
psia.  The  relative  humidity  may  now  be  calculated 


0  = 


M  5.289  x  10 


-2 


vb 


.19883 


=  .266 
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or  =  26.6% 

From  steam  tables  a  dew  point  temperature  t<}  -  21°F  may  be  interpolated 

for  pv.  =  .05289  psia.  The  density  of  the  m?xture  is  the  density  of  the 

constituents  p  =  p  +  p  and  the  partial  pressure  for  dry  air  in  state  1 
.  m  a  v 

is 

p  =  p  -  pv  =  14.78  -  .05289 
al  m  V1 

=  14.73  psia 
and  the  density  is 


Pal  (14 . 73)  (144)  ..  -2  .  3 

Pal  '  R  T_  '  '(53.T)T513)'^  "  7‘756  x  10  lb/ft 
a  a 

Pvi 

Since  0  can  also  be  written  as  0  -  — —  ,  the  density  of  the  vapour  at  1 
for  v  .  is  Pvb 

gl 

Ovi  •  »  ■  1-7”  *  10'4  lb/f'3 

and  hence  the  density  of  the  mixture  can  be  found 


p  =  p  .  +  p  =  7.756  x  10-2  +  1.733  x  lO-4 
m  al  vl 

=  7.773  x  10-2  lb/ft3 


From  Table  C2  for  test  RC30  the  volumetric  flow  rate  (V  )  of 
fresh  mixture  (m)  was 

V  =  2670  ft3-hr_1 
m 

The  mass  flow  rate  of  mixture  can  be  written  as 


m 

m 


p  V  =  2670  ft3-hr  1  x  7.773  x  10~2  lb  ft'3  x 
m  m 

5.76  x  10~2  lb  sec"1 

-2  -1  454  e 

5.76  x  10  lb  sec  x  v 

i  lb 

2.62  x  10  g’sec  1 


1  hr  _ 
3600  sec 
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This  value  compared  favourably  with  the  mass  flow  rate  readout  of  Table  C2 
26.18  g-sec~l.  The  above  calculation  is  the  only  means  of  checking  the 
mass  flow  rate  meter's  readout  which  has  been  obtained  by  aid  of  a  temper¬ 
ature  compensated  thermister  at  the  meter's  turbine  throat. 

The  mass  flow  rate  of  the  mixture  may  now  be  used  to  calculate  the 
overall  volumetric  efficiency  in  equation  (8).  For  the  calculations  of 
this  report  the  displacement  volume  (Vj)  was  not  multiplied  bv  two  as  was 
suggested  in  Reference  [2]. 


nv  ~  NV.p” 
d  m 

m  =  2.62  x  10  g‘sec  3  N  =  4000  RPM 

m 

or  (5.76  x  10"2  lb-sec-3)  p  =  7.773  lb/ft3 

m 

V.  =  18.49  in3 
d 


-2  -1 

5 . 76  x  10  lb  sec _ x  60  sec _ 

v  , „„„rev  ,n  .  3  lb  1  ftT 

4000 — —  x  18.49  m  x  7.773 — r-  x - 5- 

min  ft3  1728  in3 


=  103.8% 

The  fuel/air  ratio  of  equation  (9)  may  now  be  calculated 


F  =  ---  ,  mc  =  1.64  g-sec  1 
m  r 
a 

m  -  p  Q  ,  where  Q  -  Q 
a  a  a  a  m 

,  n~2  lb  2670  ft3  1  hr 
ma  =  7.756  x  10  -3  x  - - x  3^-— 

=  5.7  52  x  10~2  x 

sec  1  lb 

=  2. 612  x  10  g-sec  3 


F  =  - ii^i-SlSec -  =  6.279  x  10_2 

2 . 612  x  10  g-sec 

or  the  air/fuel  ratio  of  equation  (10)  is 


F  =  15. 93 
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and  the  equivalence  ratio  of  equation  (11)  is 


.  _  F(actual) 
F(stoich) 


6.279  x  10 
6.68  x  10”2 


=  .9399  . 
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A  Fuel  Engine  Evaluation  System  (FEES)  has  been  designed  to  evaluate 
the  effects  of  changes  in  fuel  composition  on  engine  performance  and 
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